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Confidence  Interval,  exponential  parameter 
censoring 


A life  test  on  Items  assumed  to  have  an  exponential  lifetime  Is  often 
designed  with  type  I censoring,  type  II  censoring  or  both.  In  type  I censoring 
n Items  are  placed  on  test  and  observed  for  a fixed  time  t*,  while  In 
type  II  censoring  the  test  terminates  with  the  failure,  where  r is 
a preassigned  Integer.  If  these  schemes  ace  combined,  the  test  terminates 
at  m1n(t*,  T-)  where  .t,,  Is  the  time,  of  the  rth  failure.  Procedures-*-*;- 


■-uiiOiAnzi.,, 


(ABSTRACT,  Continued) 

for  estimating  the  exponential  parameter  from  this  combined  scheme  were  first 
considered  by  Epstein.  Epstein  established  one  sided  confidence  Intervals  for 
this  parameter.  This  report  reviews  Epstein's  work  and  establishes  two  sided 
confidence  Intervals.  The  confidence  Intervals  are  expressed  In  terms  of  time 
on  test  and  chi-square  percentiles  with  the  degrees  of  freedom  depending  on  the 
number  of  observed  failures.  An  expression  for  the  expected  length  of  the 
confidence  Intervals  is  also  derived. 
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Reliability  acceptance  tests  are  conducted  under  a variety 
of  test  schemes  both  in  the  military  and  in  Industry.  Various 
combinations  of  time  and  failure  censoring  are  commonly  utilized. 
In  time  censoring,  n items  are  placed  on  test  and  observed 
for  a pre-chosen  time  t*  . In  failure  censoring,  n items  are 
observed  until  the  occurrence  of  the  rth  failure.  These 


procedures  are  sometimes  denoted  type  I and  type  II  censoring 
respectively.  These  two  schemes  may  be  combined  so  that  the  test 
terminates  at  mln(t*,T_)  , where  is  the  time  of  the  rth 


failure  in  the  sasqple.  Sequential  tests  are  also  used  where 
acceptance,  rejection,  and  continuation  regions  are  defined  in 
the  time-failure  plane  according  to  the  theory  of  Wald  [10].  The 
sequential  test  may  be  further  modified  by  time  censoring,  failure 
censoring,  or  both.  If  a sequential  test  is  censored  at  r^ 
failures  and  at  time  Tq  , the  resulting  test  is  known  as  a 
truncated  sequential  test  scheme.  Further  variations  are  possible 
by  monitoring  total  accumulated  time  on  test  rather  than  actual 
elapsed  time  in  any  of  the  above  schemes. 

A reliability  test  is  concerned  with  a parameter  6 in 
the  underlying  distribution  of  lifetimes.  In  addition  to  an 
accept  or  reject  decision,  it  might  be  of  interest  to  form  an 
interval  estimate  of  6 at  the  time  of  decision. 


The  objective  of  this  report  is  precisely  that.  If  the 
reliability  test  terminates  at  time  t.  , we  wish  to  form  a 
100 (1-a)  percent  confidence  interval  for  6 at  that  time. 
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I 


In  this  report  we  restrict  our  work  to  an  underlying  distribution 
of  lifetimes  which  is  exponential  with  parameter  6 . This  report 
was  motivated  by  a communication  from  the  Office  of  Naval 
Research  concerning  proposed  revisions  in  MIL-STD-781.  This 
standard  gives  detailed  reliability  acceptance  procedures  for 
military  use,  en^loying  schemes  which  combine  types  I and  II 
censoring  as  well  as  truncated  sequential  tests.  The  underlying 
distribution  is  assumed  to  be  exponential  with  parauneter  6 / 
which  is  the  mean  time  between  failures  (MTBP) . The  commtinication 
indicated  a desire  to  form  confidence  intervals  for  6 directly 
from  the  testing  procedures  in  NIL-STD-781. 


II..  Literature  Review  and  Preliminaries 

Much  of  the  existing  theory  on  reliability  testing  as  well 
as  estimation  procedures » in  the  exponential  case  was  developed 
bv  Beniamin  Epstein  and  Milton  Sobel  in  a series  of  papers  and  tech- 
nical reports.  Epstein  eventually  unified  these  papers  in  a mono- 
graph (4)  in  1960. 

In  what  follows  it  is  assumed  that  all  test  items  are 
drawn  at  random  and  independently  from  a population  where  the 

T / of  an  item  is  exponentially  distributed  with  density 

f(r/$)  ■ #Ta0,e>0  (i) 

0 f elsewhere. 


Hence  6 is  the  expected  lifetime  of  any  item  on  test.  Suppose 
n items  are  drawn  at  random  from  distribution  (1)  and  placed  on 
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life  test.  Define  to  be  the  observed  failure  time. 

If  Tin  Is  preassigned  for  type  II  censorlngr  we  have 

. Palled  items  may  not  be  replaced  (sampling 
without  replacement)  or  they  may  be  Immediately  replaced  by  a 
new  Item  from  the  same  distribution  (8aii%>ling  with  replacement) 
Under  these  schemes  the  following  results  are  well  established 


1)  In  the  nonreplacement  case  the  maximum  likelihood 


estimator  of  6 is  § « T^r  , where 


r r r 

P_  - S T.  + (n-r)T_  ■ Z (n-i+1) (T.  - T.  .)  - Z u.  . 

^ i-1  ^ ^ i-1  111  1 


Here  Tq  = 0 and  u^^  5 (n-i+1)  (t^  - is  in  fact  the 

total  accumulated  time  on  test. 

ii)  The  u^  0 i «!»...» r , are  independent  and  identi** 
cally  distributed  with  common  probability  density  function  (1). 

iii)  2Uj|^/6  is  distributed  as  a chi>square  random  variable 
with  2 degrees  of  freedom.  Consequently  2T^e  is  distributed 
as  chi**squ2ure  with  2r  degrees  of  freedom  (X^^)  • 

iv)  In  the  replacement  case  the  maximum  likelihood 

A 

estimator  of  6 is  6 ■ T^r  # where  now  ■ nt^  . 

V)  Also  in  the  replacement  case  2re/9  • 2T^e  - . 


Frosi  these  results  we  can  obtain  one  and  two  sided 


confidence  Intervals  for  6 « with  or  without  replacement.  From 
(ill)  and  (v)  we  find  that  a 100 (1-a)  percent  one  sided  confidence 
Interval  for  6 is 
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I 

i 

! 

i 


(2) 


2 2 2 

%diere  ouch  that  ® corresponding 

two-sided  interval  is 


In  type  I censoring^  where  testing  continues  until  a 
preassigned  time  t*  , confidence  intervals  are  not  as  easily 
found.  In  this  case,  with  n Itons  on  test  with  replacement, 
we  have  a Poisson  process  with  parameter  n(^)  . Hence,  if  k 
is  the  number  of  failures  observed  in  (0,t*)  , 


P(k  s r|e) 


By  a %#ell  kno«m  relationship  between  the  Poisson  and  chi-square 
distributions  we  can  write 

P(ksr|0)  • PCXj^^j  > 2nt*/0)  . 

Clearly,  if  0 s 2ntVX^  then  P(ksr)  so  . So,  if  we  do 

in  fact  observe  k«r  we  are  at  least  100 (l-o)  percent  confident 
that  2nt*/X^^2r+2  ^ ® * Likewise,  by  considering  P(k2r|0)  , 

we  may  be  100 (l-o)  percent  confident  that,  given  we  observed  k«r  , 
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As  a consequence#  one  sided  and  two  sided 


100 (1-a)  percent  confidence  intervals  for  6 , under  this  seunpling 


This  approach  to  confidence  intervals  for  the  Poisson  parameter 
was  originally  discussed  by  Garwood  [8]*  Note  that  if  r*0  , 
only  the  one*-sided  interval  is  used  since  we  have  no  information 
on  which  to  base  an  upper  limit. 

If  the  type  I scheme  proceeds  without  replacement# 
corresponding  results  are  not  easily  obtained.  In  this  case 
the  distribution  of  the  maximum  likelihood  estimator  was  shown 


by  Bartholomew  [2]  to  be  a weighted  sum  of  chi-square  Integrals. 
A closed  form  expression  for  confidence  intervals  from  this  is 
unlikely.  Epstein  used  a non-parametrlc  approach  [4]  to  this 
problem  by  considering  the  number  of  failures  observed  in  (0#t*) 
but  ignoring  the  actual  times'  of  failure.  Ne  may  improve  our 
position  considerably  here  if  we  choose  to  monitor  total  time 
on  test  rather  than  actual  elapsed  time.  If  at  some  time  t # 


k failures  have  been  observed#  the  total  time  on  test  is  seen  to  be 


£ ^4  + (n-k)t  . 

i=0  ^ 


Instead  of  an  actual  time  of  truncation  t*,  we  preassign  a total 
time  truncation  point  t'  . Again,  Epstein  [4]  has  considered 
this  approach  and  many  resultant  properties.  Under  this  scheme 
it  can  be  shown  that  we  are  observing  a Poisson  process  with 
parameter  ^ ^ length  of  actual  time  t"  . Therefore, 

following  the  seune  reasoning  that  lead  to  (4)  and  (5),  one-sided 
and  two-sided  100 (1-a)  percent  confidence  intervals  for  6 are 


"o,2r+2 


^2r+2 


^l-|r2r 


With  this  scheme,  the  objective  of  terminating  the  test  before 

an  excessive  amount  of  time  has  elapsed  may  still  be  achieved. 

A more  ccmqplex  situation  arises  when  types  I and  II  are 

combined.  In  a combined  scheme,  the  test  may  terminate  either 

with  tQ  failures  or  at  time  t*  , whichever  occurs  first.  If 

we  are  testing  Hq:  9 * 6q  against  6 ■*  6^  where  ^ ^0  * 

we  would  accept  if  min{T_  it*)  ■ t*  and  reject  if  mln(T  ,t*) 

T • It  may  be  that  r^  and  t*  are  prechosen  to  satisfy 
0 ” 

• 

a,0  level  requirements. 


Figure  1.  Possible  sample  paths  in  the  combined  test  scheme 


Additional  complexity  in  the  stopping  rule  complicates 
the  formulation  of  a confidence  interval  at  the  time  the  test 
terminates.  The  natural  Inclination  here  is  to  combine  the 
confidence  intervals  from  the  types  I and  II.  Epstein  has  done 
this  and  showed  [4]  that  the  following  rule  gives,  at  least, 

100 (1-a)  percent  one-sided  confidence  intervals  when  sampling 
with  replacement: 

'0-1  ' 

\*o,2k+2  / 
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If  we  sample  without  replacement  we  face  the  same  problems  as 
above  when  monitoring  actual  elapsed  time.  If,  instead,  we 


monitor  total  time  on  test  and  use  total  time  truncation  t'  , 
the  following  one  sided  100 (1-a)  percent  confidence  intervals 
for  6 were  also  given  by  Epstein  [4]: 


if  k«0,l. 


(9) 


Here  k is  the  number  of  failures  observed  when  total  accumu- 
lated test  time  is  t'  , and 

1^0 

t_  ■ Z + (n-r-)T  * total  time  on  teat  at  failure. 

*^0  i-1  ^ ^0  ^ 

Two  sided  intervals  are  considered  in  the  next  section. 


r 

III.  Extension  to  Two  Sided  Intervals 

Ilia.  Testing  with  Replacement 

It  would  seem  that  an  extension  from  one  sided  to  two 
sided  intervals  should  now  be  simple  and  direct.  In  fact  the 
proof  is  not  a simple  extension  from  the  one  sided  case. 

Theorem  1 gives  two  sided  intervals. 

Theorem  1.  Let  t be  the  lifetime  of  an  item  having 
distribution  f(T)  «=  (l/0)e“^'^®  , TiO  . If  n items  from 
this  distribution  are  placed  on  test  with  replacement,  and  if 
the  test  terminates  at  roin(t*,T  ) , then  the  following  rule 

gives  a two  sided  100 (1-a)  percent  confidence  interval  for  6 : 
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he  used  X*  ^ in  the  denominator  when  k » 0 . He  was  unable  to 
provide  a proof  of  his  conjecture.  To  prove  the  intervals  (1) 
are  valid  we  need  to  show  that,  wherever  6 be  in  the 

probability  is  at  least  1-a  that  6 is  covered  by  one  of  the 
intervals  (10).  The  proof  is  somewhat  tedious  and  requires 
attention  to  figure  2.  The  endpoints  of  the  intervals  for 
k * 0,1, . . . ,rQ-l  form  a partition  of  the  parameter  space.  The 
proof  consists  of  establishing  the  required  probability  of  l-*a 
for  each  sxibset  of  the  partition.  Note  first  that  the  lower 
endpoints  of  the  intervals  when  k » 0, 1,2, . . . ,rQ-l  are  ordered 
among  themselves,  the  smallest  when  k > r^-l  , the  largest  when 
k « 0 . This  is  a result  of  the  fact  that  X?  ^ . for 

any  fixed  a , emd  any  integer  k . Likewise  the  upper  endpoints 
when  k * 1, 2, . . . , rQ-1  are  also  ordered,  the  smallest  when 
k ■ Tq-I  , the  largest  when  k « 1 . The  endpoints  of  the  interval 

when  k * r,.  are  random  since  they  are  functions  of  t , the 

time  of  the  failure.  However,  < t*  implies  this  random 

lower  endpoint  is  bounded  above  by  the  fixed  lower  endpoint  of 
the  interval  corresponding  to  k « » specifically 

2 2 

2nT^  /X^y2^2r  ^ 2r  * Also  this  random  upper  endpoint 

is  bounded  above  by  the  fixed  upper  endpoint  when  X » r^  , i.e., 

2nT  /X?  -/o  I-  ^ 2nt*/X?  • Note  that  if  the  smallest 

Tq  i-a/2,2rQ  l-a/2,2rQ 


upper  endpoint  (k"rf.-l)  was  always  greater  than  the  largest 


lower  endpoint  (k  • 0) , there  would  then  be  a convenient  overall 
ordering.  Unfortunately  this  is  not  necessarily  the  case.  The 
partitioning  of  the  parameter  space  may  be  complicated  by  what 


we  shall  refer  to  as  a "cross-over".  A cross-over  occurs  when 
an  upper  limit  for  one  fixed  interval  crosses  below  the  lower 

limit  of  another  fixed  interval.  The  first  such  cross-over  will 

2 2 
occur  when  2nt*/X,  o o becomes  less  than  2nt*/X„  # 

^rQ-.^  y,  2 


that  is  when 


2x^-2  > Xa 


,2  • 


This  can  occur  if 


is  large 


enough  or,  more  likely,  if  r^  is  large  enough  for  a fixed  a 
level.  Clearly  a vast  number  of  possible  configurations  exist 
where  upper  endpoints  are  intermixed  with  lower  endpoints.  One 
such  configuration  is  illustrated  in  figure  3 below.  In  this 
illustration  a».10  amd  r^ « 13  . 
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Figure  3.  Exan^le  of  Interval  endpoints  Intermixed,  a <■.10,  rQ«13. 

All  of  these  possible  configurations  may  be  reduced  to  the  consid- 
eration of  a relatively  few  cases  as  Indicated  in  figure  2.  Clearly 
we  cannot  consider,  individually,  all  possible  subsets  of  every 
possible  configuration  of  endpoints  in  the  partition.  In  each 


conceivable  case,  however,  the  proof  would  utilize  only  a subset 


We  first  prove  two  short  lenunas  that  will  be  used  in  the 


Lemma  1.  If  xzO,  then  f(x) 


which  is  a telescoping  sum.  Hence  f'(x) 


lifetime  with  parameter  6 . That  is  f(T.)  « (l/e)e 


T.  a 0 . Then  Y » 2nT  /8  has  a chi-square  distribution  with 


To  prove  this  we  write  t 


Since  we  are  observing  a Poisson  process  with  parameter  A « n/e 

it  follows  that  T_  is  the  sum  of  independent.  Identically 

rg  Q 

distributed  random  variables  with  connon  exponential  distribution 


X z 0 . Consequently  t has  a gamma 


distribution  with  parameters  a • r^  , 0 > e/n  . If  we  let 

y ■ 2nT  /6  , the  distribution  of  Y is  easily  found  to  be  X 
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Proof  of  Theorem  1.  Case  (1)  corresponds  to  region  1 
in  figure  2. 

Let  7 represent  the  family  of  intervals  (10) » and  k 
be  the  number  of  failures  observed  in  (0,t*). 


Case  (1) : 8 2 


2nt* 
’‘i-a  2 

1 j,2 


Under  this  restriction  on  6,  P(8c7)>  P(no  failures  before  t*) 

nt* 

> e ^ , since  we  are  observing  a Poisson  process  with  parameter 

n 9 

^ . By  the  case  (1)  assunqptlon  nt*/0  i which  gives 


nt* 


-x:  o 2/2 

2 e . 


Mow  we  can  make  use  of  the  %rell  known  relationship 


P(X 


Jo  - 


(U) 


Q /2 

Thus,  with  r-0,  P(X^  > xJ_o  j)  - e ^"7'^  .1-1. 


So 


we  have  P(6c7)  2 ^**0  which  was  to  be  proved  for  case  (1) 


Case  (2) t 


i e < 


*1-5#  2h+2 


^l“f»2h 


Th  1 h 

(0  € I)  ■ ^ failures  before  t*)  ■ [ P(r  failures  in  t*) 


r«0 


nt* 


h . » (2|i)' 
•I  ® 


r-0 


rl 
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By  case  (2)  assumption#  nt*/9  ^ ^l-®,2h+2'^^  which  implies 


nt* 

h ” 6 /nt*t  r 
e (~g-) 


v2 
-X,  a 


>/2 


rl 


h « ''l-7,2h+2'‘.5,2  ^ 

^ e 2 <^l-|,2h+2/2^ 


TT 


- 1 - 


r-0 


The  inequality  follows  from  Lemma  1 and  the  equality  from  (11) 
Hence#  again#  P(ecl)  > 1-a  . 


2nt*  ^ . 2nt* 

Case  ( 3 ) : — s — s 0 < — = 


r2 


^l-|,2m 


This  is  included  in  the  more  general  case  (6)  with  d>0  and 

s > m . 


c*M  (4)1  — s e < 

*a  j 


Again#  this  is  included  in  case  (7)  with  d«0#  s>m 


Case  (5)j  < 9 < -s--  - 


^2s'*-4 


f,2s+2 


Note  here  that  we  necessarily  have  s a d-fl  . Otherwise#  if 


s - d # we  would  have  or  X?  o , , > xf 

X?  a ..  ? f,2d 

S#2d  ^ 
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which  is  clearly  not  possible.  Now  within  case  {5), 

r«-i 


p(e 


s-1 

) - I 

J r«2+] 


failures  in  t*) | « I P(r  failures  in  t*) 

r«2+l 


s-1 

Y . » (S|i)Vr. 

r-2+1  ® 


s-1 


nt* 


r 6 .nt*>r,  , 
L ® (“5“)  /r! 


r-0 


nt* 

I e ® 

r-0  ” 


X?  o 

By  case  (5)  assumption  0 > a ■— — and  so  ^ 

X*  ^ 9 2 

2« 


y2 

S-1  "*l-f,2s/2  , 

Consequently  the  first  sun  is  > 2 * 


r-0 


1 -SL 

1 2 ' 


and  P(e  € I)  > 1 - » - 


nt* 

z .~®"(S|i,Vr. 

r-0  ® 


But 


nt*/0  > X*  y,  and  so 

f,  22+2/2 


nt* 


I < I /j,r 


r-0 


r-0 


«22+2 


/rl  - ^ 


Finally  P(0€l)  >i-|-|-i-o 
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a ’^l-|,2s+2/^ 

I ® 

r-0 


^01  j , <^/2 


d 2d+2'  , 

- I e <^a 

r-O  ^,2d-»-2 


/2 


)Vrl 


Again  the  inequality  follows  from  the  restriction  on  6 which 
lollies  nt*/e  s ^1-^,  2s+2'^^  nt*/e  > 2d+2^^  * 

again  P(ec7)-  l-^-j«l-a.  It  is  noteworthy  here  to 

consider  this  case  when  d>0r  s»m»  which  is  merely  case  (4) 
as  previously  cited.  The  proof  here  relies  on  the  fact  that 
2 

6 < 2nt*/X„  . But  in  Epstein's  conjecttired  intervals^  this 

partition  point  would  be  2nt*/X^  > r which  gives  P(8  c I) 

<Xt  i 

l-j-a«l-^.  Problems  arise  if  2nt*/X^^2,2  ^ ^ 2nt*/X^  2 • 

For  exanqpler  it  can  be  shown  that  for  the  sinqple  case  with  r^  • 2 , 
with  6 in  this  small  intervalf  P(6cl)  < l-a  . 

Case  (8):  - s 0 < . 

^l-|,2v+2  ^l-|,2v 

V 


P(6  c I)  > \ P(r  failures  before  time  t*)  • 

r"d+l 
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„ nf 

I . 


r-0 


nt* 

0 

6 


I e ” (2|^)Vrl  > 
r«0 


> I • <>‘U,2V2/2)-^/rl 

r*o  2 


d "^o/2,2d+2/^ 

I e 

r-0 


(X;  /2)Vr! 

|,2d+2 


2 2 

Again  this  follows  from  nt*/0  s ^i-®,2v+2  nt*/0  > X^^2,2d+2 


So,  for  case  (8),  P(0el)  > 1 - 1-a  . 


I' 

0 

I " 2 


Case  (9) : 


2nt* 


|/2w+2 


i e < 


2nt* 


v2  « 

^l-|,2rQ-2 


Mote  that  now  we  must  include  the  possibility  of  coverage  by  the 
random  interval  corresponding  to  h - r^  . 


^0-1 


2nT. 


P(0  € I)  » I P(r  failures  before  t*)  + ^ 


r«w 


Xl«“ 

1-y/  2rQ 


^0*1  -2^ 

I e ® (2|^-)Vrl  + P( 


0X2  o 


r-w 


1-f , 2rQ 

Tn 


STS  t*)  . 

*^0 


It  follows  from  Lemma  2 that  we  can  write  this  second  term  as 


P( « T,  s t*)  - P(xJ_o  ,,  s Y s 2*^) 


F„<2ntV0)  - # - 


Tq-I  nt* 

l-f-P(x’r^>2!^)  - 1-f-  I •^(S|^>Vrl 


and  so  nt*/e  > X 


0 < 2nt*/X 


which  then  implies  P(e  e I)  > 1 


2nt* 


By  the  same  reasoning  as  in  case  (9) » 

V 

• 1 a ,nt*t  r 


Case  (11) : 6 < 


t»2ro 


p(e  € I) 


2nT. 


2nT. 


$ e ^ 


l'2ro 


^1——  2r 


ex 


l-a/2,2r, 


ex 


a/2,2r. 


2n 


s T. 


2n 


^^^1— — 2r 
■*■  2'^‘0 


i X 


2r, 


i X 


1 - a 


}>2x, 


These  cases  cover  any  possible  subset  of  every  partition.  A formal 
proof  of  this  claim  is  found  in  Appendix  1.  The  proof  of  Theorem  1 
is  complete. 

mb.  Expected  Interval  Length 

It  is  of  some  interest  now  to  derive  an  expression  for 
the  expected  length  of  the  confidence  interval  under  this  scheme. 
There  is  a positive  probability  of  observing  no  failures,  in  which 
case  the  interval  length  is  Infinite.  To  avoid  this  situation 
we  will  condition  on  one  or  more  failures.  Denote  the  length  of 
the  interval  by  L and  consider  the  following  two  possibilltlest 
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If  T > t*  then  L = 2nt*  -5-= 


2k+2_ 


, say. 


Here  k is  the  number  of  failures  observed  with  1 ^ k 


If  T it*  then  L 


y2 

2r 

2'^^0 


[ » CT  , say. 

^0  ^0 


Note  that  in  the  first  case  the  length  is  discrete  while  in  the 
second  it  is  a continuous  remdom  variable.  Consequently  the 
unconditional  distribution  of  L is  of  the  mixed  type  with  jumps 
at  ^k^  ^*^f2f***^^Q*^  • So# 


P(i)  « / f(t)dt  + I P(k  failures  in  t*)  if  A < ct* 
0 kcA 


/ f(t)dt  + I P(k  failures  in  t*)  if  A i ct*  , 
0 kcA 


where  f is  the  density  function  of  t and  A >■  {k:  A.  s A}  . 

^0  * 

From  Lemma  2 we  have  that  t ~ gamma  (a^r.#  B * Q/n)  # and 

_nt* 

P(k  failures)  ■ e ® . We  want  E(L|L<»)  which  is 

E(L|k7<0)  since  L«a>  if  and  only  if  k"0  . We  need  now  the 
conditional  distribution  of  L . We  can  write 


P(L^A|L<«>)  * P(L  s A n L < ••)/P(k  ^0)  . Now  if  A is  finite# 
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(LSI)  c (L  < «)  inrhich  then  implies  P(Ls£  n l < *)  » P(Lsl)  . 

On  the  other  hand,  if  I is  infinite,  (L < •)  c (L  s t)  and 
P(Lsl  n L<<»)  « P(L<“)  . Thus 

F(£|l<«)  « finite 


» 1 


if  I is  infinite. 


Consequently,  E(L|)c,<0)  - / idF{i|kyO)  - _^^/*ldP(£)  , 


nt* 

where  POcf^O)  « 1-e  ® . Thus, 


E(Llkf<0) 


P()tj<0) 


ct* 


I fail  in  t*)  + / £g(£)d£ 

k-1  0 


1 

P{k><0) 


r 

I KP(k  fail  in  t*)  + / CT^  f(T^  )dT 


0 *0  *"0 


To  evaluate  the  integral* 


t* 

c / yf(y)dy 
0 


- c / 


t* 


.-f 


0 r(r,)(2)  0 


rUjI  (J)"® 


-,t* 


ro. 

* i ^ V+l 

k-0  <rQ-k)l(-|)*^^ 


y«0 
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n 


1 - 


*"0 

I 

k-0 


nt* 


(Htl  k 
' 0 ’ 


kl 


Thus, 


E{L|kf<0) 


nt* 


1-e 


I 

k«l 


nt* 

0 


(— )^ 
' 0 * 


kl 


0 c r. 


n 


k-tQ+l 


e ' 0 ' 
ki 


where 


2nt* 


v2 

2k 


^2k+2 


Several  questions  may  come  to  mind  concerning  the  expression  for 
E(L|k^0)  . For  instance,  if  we  take  the  limit  as  r^  * in  (12), 
do  we  obtain  the  E(L|kf<0)  for  the  type  I censoring  scheme? 

If  we  take  the  limit  as  t*  , do  we  obtain  the  E(L)  for 
the  type  II  censoring  scheme?  It  can  be  shoim  quite  easily  that 
the  answer  to  both  questions  is  yes.  In  the  preceding  derivations 
we  have  placed  no  restrictions  on  t*  and  r^  other  than  to  say 
they  are  censoring,  or  truncation,  values.  In  practice  the  values 
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of  t*  and  are  determined  by  the  test  procedture  which 

provides  the  data  for  estimation.  MIL-STD-781  is  concerned  with 

testing  HqI  0 * Qq  against  6 * where  ^ ®0  * • 

test  combining  type  I and  type  II  censoring  with  P(reject  «a 

2 

is  desired,  Epstein  {5]  has  sho%m  that  t*  * 6„X.  ..  /2n  . 

u A—a# 

If  in  addition  we  require  P (accept  s 8 , Epstein  has 

shown  this  may  be  achieved  if  we  choose  r^  to  be  the  smallest 


integer  such  that 


e,2r, 


l-a,2r, 


Then  we  may  decrease  8 if  we  increase  r^  . It  seems  reasonable 
that  an  increase  in  r^  , or  a decrease  in  8 » would  result  in 
a corresponding  decrease  in  E(l.|]c^0)  . An  attempt  at  a direct 
analytic  proof  using  (12)  was  not  successful.  We  then  programmed 
(12)  and  generated  E(L|k7<0)  for  various  combinations  of  6,0^ 
and  a . The  value  of  rQ  was  increased  for  each  combination. 

The  results,  somewhat  surprisingly,  show  that  for  some  combinations 


of  0r0A  And  a , the  expected  length  incr< 


initially  before 


it  does  decrease  with  increasing  r^  . Partial  results  of  this 
computer  work  are  given  in  Table  1.  We  have,  thus  far,  been 
unable  to  provide  a satisfactory  explanation  of  this  behavior. 


r 

I 


1.00 


e. 


3.00 


a ■ 0.05 


0.10 


• 

2 

E(L) 

m 

25.7057 

R 

m 

2 

E(L) 

m 

14.6718 

m 

3 

E(L) 

m 

26.1465 

R 

m 

3 

E(L) 

m 

11.1729 

m 

4 

B(L) 

m 

15.0962 

R 

m 

4 

B(L) 

m 

5.6256 

m 

5 

B(L) 

m 

6.9176 

R 

m 

5 

B(L) 

m 

2.8950 

R 

m 

6 

E(L) 

m 

3.4234 

R 

m 

6 

E(L) 

m 

1.9311 

m 

7 

E(L) 

m 

2.2406 

R 

m 

7 

BCD 

m 

1.5818 

m 

8 

E(L) 

m 

1.8238 

R 

m 

8 

E(L) 

m 

1.4094 

m 

9 

E(L) 

m 

1.6299 

R 

m 

9 

Ea) 

m 

1.2954 

m 

10 

E(L) 

m 

1.5034 

R 

m 

10 

E(L) 

m 

1.2071 

m 

11 

B(L) 

m 

1.4068 

R 

m 

11 

E(L) 

m 

1.1351 

m 

12 

E(L) 

m 

1. 3266 

R 

m 

12 

E(L) 

m 

1.0745 

R 

■ 

13 

E(L) 

m 

1.2590 

R 

m 

13 

E(L) 

m 

1.0223 

a 

14 

E(L) 

1.2001 

R 

m 

14 

E(L) 

m 

0.9771 

m 

15 

BCD 

m 

1.1490 

R 

m 

15 

E(L) 

m 

0.9374 

6 - 2.00 

0-  » 3.00 

a 

• 0.05 

a 

■ 0.10 

R 

2 

E(L) 

m 

33.6920 

R 

2 

E(L) 

m 

22.2509 

R 

3 

E(L) 

m 

53.1913 

R 

3 

E(L) 

m 

28.8162 

R 

4 

E(L) 

m 

54.6214 

R 

4 

E(L) 

m 

25.6829 

R 

5 

E(L) 

m 

44.3260 

R 

5 

E(L) 

m 

18.7492 

R 

6 

B(L) 

m 

30.9006 

R 

6 

E(L) 

m 

12.2806 

R 

7 

E(L) 

m 

19.6450 

R 

7 

ECL) 

m 

7.7827 

R 

8 

E(L) 

m 

12.0050 

R 

8 

B(L) 

m 

5.1069 

R 

9 

E(L) 

m 

7.4825 

R 

9 

E(L) 

m 

3.6492 

R 

10 

E(L) 

m 

5.0097 

R 

10 

E(L) 

m 

2.8815 

R 

11 

E(L) 

m 

3.7223 

R 

U 

E(L) 

m 

2.4700 

R 

12 

E(L) 

m 

3.0517 

R 

12 

B(L) 

m 

2.2326 

R 

13 

E(L) 

m 

2.6882 

R 

13 

E(L) 

m 

2.0790 

R 

14 

E(L) 

m 

2.4717 

R 

14 

Ea) 

m 

1.9683 

R 

15 

B(L) 

m 

2.3277 

R 

15 

sa) 

m 

1.8805 

e - 4.00 


3.00 


a * 0.05 


a • 0.10 


2 

Ea) 

- 38.6661 

R - 

2 

sa)  - 

3 

Ea) 

• 73.4801 

R - 

3 

Ea)  - 

4 

Ea) 

- 97.9833 

R - 

4 

Ea)  - 

5 

Ea) 

-109.0797 

5 

Ea)  - 

6 

Ea) 

-108.3810 

R - 

6 

Ea)  - 

7 

Ea) 

- 99.3520 

R - 

7 

Ea)  - 

8 

Ea) 

- 85.6724 

R - 

8 

la)  - 

9 

Ea) 

- 70.4874 

R - 

9 

KL)  - 

10 

Ea) 

- 55.8853 

R - 

10 

Ea)  - 

11 

sa) 

- 43.1227 

11 

sa)  - 

12 

Ea) 

- 32.6525 

R - 

12 

la)  - 

13 

la) 

- 24.4603 

R - 

13 

■a)  - 

14 

sa) 

- 18.3067 

R - 

14 

Ba)  - 

15 

sa) 

- 13.8143 

R - 

15 

Ba)  - 

27.0906 

44.5701 

54.0796 

56.0089 

52.4240 

45.7010 

37.8212 

30.1389 

23.4156 

17.9371 

13.6968 

10.5338 

8.2472 

6.62S6 
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Table  I.  Expected  Lengths  for  Increasing  Truncation  Values, R 
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IIIc.  Testing  without  Replacement 

The  second  testing  scheme  we  are  considering  terminates 
at  ain(T  ft"),  where  failed  items  are  not  replaced.  Recall 
that  the  truncation  value  t'  now  refers  to  the  total  accumu- 
lated test  time.  At  the  time  of  decision  the  following  rule 
gives  a 100 (1-a)  percent  confidence  interval  for  6 : 


2t' 


a , 


if  k-0 


2t‘ 


2t' 

' "5 

X,_a 


' 2k-*-2 


if  k 


Tq”!# 


2T 


2T 


» 2 

x:_o 


fjL 


If  k 


^ i*  the  number  of  failures  observed  when  the  total  time 
on  test  reaches  t'  , and  T is  the  total  time  on  test  at  the 
time  of  the  r^th  failure.  This  rule  is  nearly  Identical  to  that 
conjectured  by  Epstein  [4],  with  again  the  only  change  being  the 
use  of  a/2  rather  than  a when  k-0  . The  proof  of  (13) 
becomes  simple  when  we  realise  that  this  schema  may  be  viewed  as 


(13) 


a Poisson  process  with  parameter  1/6  which  we  observe  for 


actual  time  t'  . To  see  this,  let  Tj^  be  the  total  time  on 
test  when  the  item  fails.  Then 


k k 

“ I T + (n-k)T.  - I (n-i+1) (T . - T . ,) 

K i_i  1 K 1 1-1 


where  we  have  previously  noted  the  are  independent  and 

identically  distributed  as  exponential  (j)  . It  follows  that 
Tj^  ~ gamma  (a«k,  6«6)  . We  siay  trrite 

P{k  failures  in  total  time  t')  ■ -P(Tjj>t')  . From 

direct  integration  we  find  that 

k -V  t'  k 

- I • (V)  Ai 

r"0  * 


P(T,j  > t') 


k-1 

I 


r-0 


e*^(^)*‘Al 


t' 

"T  k 

and  P(k  failures)  • e (-g-)  Al  • Hence,  to  prove  (13)  we  need 

only  duplicate  the  proof  of  (10)  substituting  t"  for  t*  and 
1/6  for  n/6  . The  derivation  of  B(L|kffO)  for  this  scheme 
follows  in  the  same  manner.  The  result  is 


E(L|kf<0) 


1-e 


I 

k-1 


TT 


k-r,.-»-l 


TT 
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Appendix  I 

We  now  give  a proof  that  the  cases  considered  in  Theorem  1 
do  in  fact  cover  every  subset  of  the  possible  partitions  of  the 
parameter  space  determined  by  the  fixed  interval  endpoints. 

With  the  exception  of  cases  (1)  and  (11)  (regions  1 and  11 
in  figure  2),  each  possible  subset  in  the  partition  will  be  bounded 
above  and  below  by  a left  or  right  interval  endpoint.  Conseqxiently 
there  are  only  four  possible  subsets t 

(i)  ^ (ii)  - ] (iii)  ^ (iv)  

We  consider  each  of  these  and  show  how  each  is  covered  in  Theorem  1. 

(1) 

In  the  most  general  situation  this  is  case  (5)  of  Theorem  1.  The 
arrows  Indicate  the  bounds  on  e which  are  necessary  for  the 
proof  of  case  (5) . It  is  possible  that  there  is  no  right  endpoint 
below  the  left  endpoint  determined  by  k ■ s-t-l  . In  this  event  we 
have  s>l  • rn**!  and  we  are  then  in  case  (10)  of  Theorem  1. 

15 

(11) 


In  general  this  is  case  (6)  of  Theorem  1.  Again  the  arrows 
indicate  the  boundary  points  necessary  for  its  proof.  If  there 
is  no  right  endpoint  corresponding  to  . k ■ s**-!  * it  must  be  that 


32 


1 and  we  have  case  (9)  with  d«w  . If  there  la  no  left 


endpoint  Indicated  by  d-1  we  have  the  situation  where  d>0 


and  8-m  which  is  case  (3) . If  neither  of  the  points  indicated 
by  84-1  and  d-1  are  found,  we  must  have  the  situation  %fhere 


no  crossover  occurs.  Then  we  are  in  case  (9)  with  W"0 
d+1  8+1  d » 

(ill)  ]— e: P—  I- 


This  is  case  (7)  in  general  and  its  proof  is  dependent  only  on 
the  actual  boundaries  of  the  subset.  No  other  considerations 


In  general  this  is  case  (8)  of  Theorem  1.  If  the  endpoint 
Indicated  by  k > s+1  is  not  present,  we  must  be  in  case  (2) 
wl  th  d+1  • 0 and  v*h  • 


By  adding  regions  (1)  and  (11)  in  figure  2,  we  see  that 
Theorem  1 covers  all  possible  subsets.  Cases  (3)  and  (4)  are 
included  only  for  clarity  in  figure  2 and  may  be  omitted,  since 
they  are  special  cases  of  (6)  and  (7) . 
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